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We use a self-consistent chiral-hydrodynamic formalism which combines the linear 
o model with second-order hydrodynamics in 2 + 1 dimensions to compute the spec- 
trum of thermal photons produced in Au+Au collisions at yjs jvjv = 200 GeV. The 
temperature-dependent shear viscosity of the model, rj, is calculated from the linearized 
Boltzmann equation. We compare the results obtained in the chiral-hydrodynamic model 
to those obtained in the second-order theory with a Lattice QCD equation of state and 
a temperature-independent value of n/s. We find that the thermal photon production is 
significantly larger in the latter model due to a slower evolution and larger dissipative 
effects. 
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1. Introduction 

In this paper we focus on the spectrum of thermal photons produced during the 
expansion of the fireball created in Au+Au collisions at ^/snn = 200 GeV. Photons 
created in the interior of the fireball pass through it with out any int eraction, thus 
giving information on the properties of bulk nuclear matter! ' 1 ^ 1 4 | o | b | Thermal pho- 
ton sp ectra have been calculated within the framework of different fluid dynamics 
models L T | 6 | 5 | 3 | 4z | ag wg vj ag a hybrid kinetic- hydrodynamic model Here we employ 
a self-consistent chiral-hydrodynamic model combining second-order (SOT) hydro- 
dynamics and the linear a model (LSM), in which the temperature-dependent r//s 
is computed from kinetic theory El. We compare the spectrum of thermal photons 
with that obtained from the same SOT with a Lattice QCD-inspired equation of 
state (EOS) and a temperature-independent value of r)/s. The latter model is com- 
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monly employed in the study of hadronic observables of heavy ion collisions. The 
purpose of this contribution is to quantify the differences in thermal photon spectra 
that arise from differences in these two fluid dynamical models when they are taken 



as a whole. In other words, no attempt is made to disentangle the impact of EOS 
or the temperature d« 
study along this line). 



or the temperature dependence of rj/s on the thermal photon spectra (see SI for a 



2. Chiral-hydrodynamic model 

We now describe the coupling between the chiral fields and the quark fluid, the 
latter described by second order boost invariant fluid dynamics (see l^l^la for further 
details). 

The classical equations of motion of the LSM are given by 

D,D^ a + ^ = -gp a (1) 

where <f> — (a, if) and 

Pa = 9<f>ad q I ff 1 =fo(p»,xn (2) 

/o(p M , a? M ) is the Fermi-Dirac function and U is the potential exhibiting chiral sym- 
metry breaking. These are solved selfconsistently together with the conservation 
equations for the combined system of chiral fields plus quark fluid 

De = -(e + p)V M u" + W v a^ + q{p s Da + p ps ■ Dtt) 

1 (3) 
(e + p)Du i = -(g^dje - u l u a d a e) - A^IF' 9 + g( Ps V l a + p ps ■ W) . 

where (p s ,p ps ) = p a , D = u^D^ is the comoving time derivative, V M = A^D 01 
is the spatial gradient, and a^ v is shear tensor. The evolution equation f or 
contains terms up to second order in velocity gradients and it is given in 12141. In 
order to estimate the temperature dependence of the shear viscosity r](T) in the 
LSM, we adapt the method described in ^ and employ the linearized Boltzmann 
equation in the relaxation time approximation (see ^j. 

To compute the spectrum of thermal photons created during the evolution of 
the fireball we consider the processes of Compton scattering, qq annihilation and 
bremsstrahlung in the QGP phase, and irir —> p-f, up — > n-f and p —> ttttj in the 
hadron phase ™ as described in 121. In all cases, the nonequilibrium distribution 
function of the quarks, /(x M ,p M ), is calculated from Grad's ansatz. 

3. Results 

We now go over to present our results. As a typical example of a conformal fluid 
we take = 2(2 — In 2)-^ and Ai = — ^ corresponding to the supersymmetric 
Super- Yang-Mills plasma. We use a 13 fm x 13 fm transverse plane, and set the 
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impact parameter to b = 7 fm. The initialization time is set to To = 1 fm/c, and we 
use u x = = and (H xx , H xv , H yv ) = as initial values. The initial energy density 
profile is obtained from Glauber's model, with a temperature Tj = 330 MeV at the 
center of the fireball. We take 7?(ro,r) = and a(T ,r) = f n [l — e~( r / r °) ] with 
7'o = 9 fm as initial values for the chiral fields. We use the isothermal Cooper-Frye 
freeze-out prescription with a freeze-out temperature Tp = 130 MeV. The critical 
temperature is set to T c — 170 MeV. 

The upper panel of Figure Q] shows the square of the speed of sound c 2 = dp/de 
corresponding to the LSM with g = 3.2, w hich yields a smooth crossover, and to 
the EOS calculated by Laine and Schroder which connects a high-order weak- 
coupling perturbative QCD calculation at high temperatures to a hadron resonance 
gas at low temperatures via an analytic crossover. Recent works have shown that the 
comparison of hydrodyna mic s imulations with RHIC data favor a smooth crossover 
over a first-order transitional^, and this is the reason for considering a crossover (and 
not a first-order phase transition) in the LSM. There is no appreciable change in the 
photon spectrum provided g ~ 3.2, but if g > 3.4 or g < 3 the dip in c 2 at T c becomes 
too sharp or fades away, respectively, resulting in a temperature dependence that 
strongly disagrees with the one obtained in Lattice QCD calculations. As it can 
be seen from Figure [l] for large temperatures the conformal limit c 2 — 1/3 is 
reached, while for low values c 2 goes to zero. Note also that the transition region is 
significantly broader in the LQCD EOS than in the LSM. 

The lower panel of Figure Q] shows rj/s as a function of temperature for the LSM 
with g — 3.2. It is seen that the ratio remains small (< 0.1) in the QGP phase and 
starts increasing abruptely for T < T c . For this reason we must impose a cut-off 
in the value of rj/s; we take rj/s < 0.4. For the SOT with the LQCD EOS we 
consider a temperature-independent value rj/s = 0.11, which corresponds to the 
averaged value throughout the evolution. 

Figure [2] shows the different contributions to the total thermal photon spectrum 
coming from the QGP and the hadronic phases, for the LSM and the SOT with the 
Lattice QCD EOS. It is seen that in all cases the photon spectrum is considerably 
smaller in the LSM, the difference being larger in the hadronic phase. The difference 
between the spectra in both models is practically independent of pt, except at 
very low values where the difference becomes slightly larger. In both models, the 
contribution of the hadronic phase to the photon spectrum is dominant over that 
coming from the QGP phase. 

There are two reasons that make the spectrum obtained from the chiral model 
smaller than the corresponding to the Lattice QCD EOS. First, the hydrodynamic 
evolution is much faster in the LSM model, essentially because the speed of sound 
is on average larger than the one corresponding to the Lattice QCD EOS. This is 

a Any hydrodynamic theory necessarily breaks down when viscosity is large enough. The value of 
nls at which it is sensible to impose this cut-off can be estimated in (n/s)\ c ~ 0.4 — 0.5 (see e.g. 
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Fig. 1. (Color online) Speed of sound of the LSM and that of the model EOS of Laine and 
Schroder 15 (upper panel); r/(T) / s(T) in the LSM (lower panel). The value of the coupling constant 
is g = 3.2. 

illustrated by the fact that the freeze-out temperature is reached in ~ 6 fm/c in the 
LSM and in ~ 10 fm/c, which represents a significant difference. The other reason 
is that the shear tensor is on average significantly smaller in the LSM, leading to a 
smaller contribution of dissipative effects to the total photon spectrum. 



4. Conclusions 

Using a simple self-consistent chiral-hydrodynamic model based on the linear a 
model and 2 + 1 conformal second-order hydrodynamics, we calculated the spec- 
trum of thermal photons including nonequilibrium corrections produced during the 
evolution of the fireball created in Au+Au ultrarelativistic collisions. Compared to 
the results obtained with a commonly used second-order fluid dynamic model with 
an equation of state inspired in Lattice QCD and a temperature independent rj/s, we 
find that the photon production is significantly smaller in the chiral-hydrodynamic 
model, due to a faster evolution and to smaller dissipative corrections to the nonequi- 
librium distribution function. 
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Fig. 2. (Color online) Contributions to the thermal photon spectra from the QGP and hadronic 
phases for the LSM with g = 3.2 and the SOT with Lattice QCD EOS. 
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